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ABSTRACT 

A numerical study of the description of motion of system of particles in different resisting media was carried out. The paper 
transforms the time dependent Newton's equation of motion in a resisting medium into a space dependent form which necessitated 
the use of classical Runge-Kutta method of order four to determine the solution of the model. The work which focused on resisting 
medium proportional to exponent of velocity, showed that increase in the mass of a particle in an air resisting medium results in a 
corresponding increase in the velocity of the particle while an increase in the density of the resisting medium led to a decrease in 
the velocity of the particle which is not different in the relativistic medium. An algorithm on the application of Runge-Kutta method 
of order four is also proposed to tackle models of this nature. 


PACS (2008): 45.10.Hj;45.20.D-;45.40.Gj;45.40. -f 
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1. INTRODUCTION 


Generally, an object or particle thrown upward or allowed to fall from a height is acted upon by its weight and a host of other forces 


as well. Essential classes of forces for our consideration are those that tend to oppose the motion of an object or particles. These 
classes of forces arise because of motion in some media especially fluids referred to as resisting, damping or dissipative forces and 
the corresponding media is said to be a resisting media. For dynamic systems, a physically valid formula expressing the resistive 
force as a function of the velocity is yet to be ascertained [1] or the resistive force can have complicated velocity dependence [2]. 
However, experiments show that for low velocities, the resisting force is in magnitude proportional to the velocity. But this simple 
relationship breaks down before we reach those velocities which are of interest in ballistics. For low ballistics, the resistive force 
varies proportionally to the square of the velocity, but the law again breaks down when the velocity of the particle or object 
approaches the velocity of sound. The law of dependence is then complicated and can only be represented by experimentally 
obtained values. Therefore, a simple relationship for particle trajectories in a resistant medium is still a challenge till date. This 
difficulty, posed to researchers to propose a law or working relationship of particles in a resistant medium have not in any way slow 
down the investigation of particles motions in a resistant medium using different approximations and assumptions [2]. A study on 
the effect of Hartmann number on the part of a particle in a resistant medium was carried out by [3] wherein the conditions of the 
resistive force proportional to both the velocity and square of the velocity were examined and far reaching conclusions made. [4], 
studied the motion of a particle in a resisting medium using fractional approach and proved that, there exist a relationship between 
the coordinate and the fractional component in a given system through a physical parameter. [5], investigated the motion of a 
particle through a resisting medium of variable density. In the study, he applied a resistance proportional to the product of 
instantaneous velocity and air density and opined that the density diminishes exponentially with the altitude. The vertical projection 
in a resisting medium: reflections on observations of Mersenne, was critically examined by [6]. The study validates empirical 
observations of Marin Mersenne on timing of vertically-launched projectiles for a general mathematical model of resistance. In a 
study, [7], proposed a fractional differential equation describing the behaviour of two dimensional projectile motions in a resisting 
medium and opined that the trajectories of the projectile in the fractional approach are always less than that of the classical case. 
Studies of [8, 10], examined the effect of magnetic field on projectile motion in a resisting medium and reported that, increase in 
magnetic field result in a decrease in the velocity of the particle. [11], included electric field to the study of [8] and deduced similar 
findings. [12], in the study of horizontal and vertical projectile motion in a resistant medium subject to varying path angles and 
speed, showed that increase in path angle decrease the projectile motion. Our aim is to compare the behaviour of system of 
particles in an air resisting medium considering the non relativistic case and the relativistic case with a view to proposing a model to 
describe the different velocity dependent resisting function. 


Classical Runge-Kutta Method of Fourth Order 
The Runge-Kutta formula involves a weighted average of values of fly,¥) taken at different points in the interval y, < y< y 


n+l 


. It is given by [13] 


= v, ca (ku a 2k 49 ae 2K + kis) (1) 
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2. FORMALISM 


A system of particles of mass (m) moving with an initial velocity 10m/s in a resisting medium whose steady resistance to the 


motion of the particle acts instantaneously in the opposite direction to that of the particles motion. The mass of the particle 
experiences, the weight downward and the resisting force upward. The study assumed that; 

The resistance medium is air. 

The acceleration due to gravity is constant over the range of motion and directed downwards. 

The rotation of the earth does not affect the motion of the particle. 

The resistance or damping to the motion is a continuous function of the velocity satisfying natural physical requirements. 

The upward buoyant force is neglected 


With these assumptions, the Newton's law of motion in the presence of a resistance proportional to the velocity is given as 


a= 
md*r 1 4 
= mg —— B, pAv (2) 
dt* 2 
where m is mass of particle, 7 is position vector, g is acceleration due to gravity, B, is a positive dimensionless proportional 


constant that specifies the strength of the resisting force, t is time, 9 is density of air, A is cross-sectional area of the object and V 


is velocity vector. 


0.5 8, Av? 
v = 10m/s 
Vv 
mg mg 


Figure 1 Model and coordinate system of the problem 


Equation (2), can be written in terms of vertical space coordinate as 


(3) 


dv 2S BPA 
dy v 2m 


with the initial condition v(O) =10m/s 


3. METHOD OF SOLUTION 


Case 1 

Non relativistic 

For comparatively low velocities as observed in the case of damped simple harmonic oscillator and parachutist, equation (3) is a 
model for these descriptions. The solution to equation (3) is difficult to tackle analytically. The classical Runge-Kutta method of order 


four as stated in equation (1) is applied with a step size (h) = 0.1. Other parameters used are B, = 0.5, m(kg) = 0.145, 0.245, 0.345, 
0.445, 0.545. g =9.81m/s*, A=0.0042m* p =1.2754kg/m*> c=3x10* m/s 
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Table 1 Table 2 

m(kg) v(m/s) 

0.145 | 10.86536900 plkg / m*) wwe) 

0.245 | 10.86813181 1.2754 10.8817754 

0.345 | 10.87008446 23754 10.88164244 

0.445 | 10.87176352 3.2754 10.88150701 

0.545 | 10.87273650 4.2754 10.88137564 
5.2754 10.88122539 


Relativistic case 
For the relativistic medium, equation (3) takes the form 


di _g_ BpAyi-# 


dy v 2m, 


(4) 


Vv 
where B = —, Cis velocity of light, 7p is rest mass 
Cc 


The result of the relativistic case is not different from the non relativistic case. 


Case 2 

Non relativistic case 

For ballistics which is concerned with projectile fired upward from a cannon or small arms and free flight of bombs as well as flight 
of rockets, the resistance or damping is proportional to the square of velocity. Equation (3) in this situation, transform into 


dv_ 8 _ PB, pAv 
dy v 2m 


Using the same parameters and method, the results is tabulated as 


(5) 


Table 3 
ko) al Table 4 

m v(m/s 

a 10.88024478 plkg 1m?) | vorvs 

0.245 10.88094341 1.2754 10.88024478 

0.345 10.88123411 2.2754 10.87827250 

0.445 10.88139407 3.2754 10.87757757 

0.545 10.88149514 4.2754 10.87610770 
5.2754 10.87490364 


Relativistic case 
For the relativistic medium, equation (3) takes the form 


dig B,pAvJi-B 


dy v 2m 


(6) 


The result of the relativistic case is not different from the non relativistic case. 
Case 3 


Non relativistic case 
For sounds and others, where the exponent of the velocity is higher than or equal to three, the model transform to 
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dig Bypav™ 
dy v 2m, 


n=3 (7) 


Relativistic case 
In the relativistic medium, the model takes the form 


dig Bpav J1-B 
dy v 2m, 


The solution of equations (7) and (8) can be implemented by the algorithm below 


(8) 


Algorithm 

Define f(¥,¥) 

Input initial values Y. and V, 

Input step size h and number of steps n 
Output Y, and Vy 

For j from 1 to n do 

k= f(3.¥) 

k, = f(¥+0.5*h,5+0.5*h*k,) 
k, = f(¥+0.5*h,¥+0.5*h*k, ) 
k, = f(v+hv+h*k,) 


Wi =F +E) 420k 420k $44) 


y=yth 
Output y and Vv 
End 


The algorithm can easily be implemented on FORTRAN 9 and Python programming languages. 


4. DISCUSSION 

Case 1 

The mass of a baseball is chosen as the starting mass and subsequently, increase in the mass corresponds to an increase in the 
velocity through the resisting medium as shown in table 1. An increase in the density of the resisting medium also led to a decrease 
in the velocity of the baseball because increase in resisting medium has been enhanced as tabulated in table 2. 


Case 2 

Table 3, showed a relationship between the mass of a particle and its velocity in a resisting medium where the resistive force is 
proportional to the square of the velocity. For a typical mass of a baseball as used, its increase result in a corresponding increase in 
the velocity of the particle. This observation remains the same in the relativistic case. The equation (equation 6) shows that as long 
as v << c, theoretically, the observation in both the relativistic case and the non relativistic case shall remain the same. Table 4 also 
showed that varying the density of the medium, decreases the velocity of the particle owing to increase resistance of the medium. 


Case 3 

The pattern of relationship that exist between mass of the baseball in a resisting medium and its velocity and that of density and 
velocity in case 1 and case 2 is not different from case 3, however, owing to the exponent n, we proposed a mathematical model and 
a numerical solution using an algorithm. This will serve as a generalized solution of this nature irrespective of the value of n. 
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5. CONCLUSION 


One of the important consequences of special relativity is that mass increases with velocity. This assertion was buttressed in this 
work. Also, 0.5 is chosen for the dimensionless empirical quantity for sphere called the drag force coefficient which baseball fit into. 
From observation and some literatures cited, it is observed that objects under which resistive force is proportional to the square of 
velocity are generally larger than objects under which resistive force is proportional to the velocity. Generally, we observed that the 
space dependent form of Newton's equation of motion of particles in a resisting medium describes better the path of a particle and 
the effect of varying other parameters or factors over a long time. Above all, no significant difference between the results of 
relativistic case and that of non relativistic case of vertical motion of a particle in a resisting medium was observed in this study. 
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